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Abstract 



We discuss the renormalization of different definitions of quark masses in the Wil- 
son and the tree-level improved SW-Clover fermionic action. Using perturbative 
and non-perturbative renormalization constants, we study quark masses in the MS 
scheme from Lattice QCD in the quenched approximation at /3 = 6.0, (3 = 6.2 and 
(3 = 6.4 for both actions. A reliable extrapolation to the continuum limit is however 
not yet possible. The most reliable results are: m MS (2GeV) = 5.7± 0.1 ± 0.8 MeV , 
mf^(2 GeV ) = 130 ± 2 ± 18 MeV and m|^(2 GeV ) = 1662 ± 30 ± 230 MeV. 
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1 Introduction 



Quark masses are among the least known fundamental parameters of the Stan- 
dard Model. Due to confinement, they cannot be measured directly and our 
knowledge of these quantities relies on techniques like Chiral perturbation the- 
ory (ChPT) [1], QCD Sum Rules (QCDSR) [2]- [5] and Lattice QCD (LQCD) 
[6]- [12]. ChPT gives rather precise determinations of ratios of quark masses 
while QCDSR and LQCD determine their absolute values. Moreover LQCD 
does not require model parameters or ad-hoc assumptions. Each technique 
suffers from different sources of errors that should be carefully studied. In the 
most recent LQCD simulations the main errors are due to the quenched ap- 
proximation (i), the reach of the continuum limit (ii) and a correct matching of 
the lattice quantities to the continuum ones (iii). In this paper we try to over- 
come as much as possible the last problem. We define quark masses through 
the Vector and Axial Vector Ward Identities and discuss the renormalization 
procedure. We compute light, strange and charm quark masses with pertur- 
bative (PT) and non-perturbative (NP) renormalization constants (RC) using 
the Wilson and the tree-level improved Sheikholeslami-Wohlert (SW) Clover 
action [13], [14] which, in the following, we will be referring to as Improved 
Action. We also try to estimate the remaining overall systematic uncertainty, 
mainly due to (i) and (ii), comparing the chiral behaviour of the pseudoscalar 
and vector meson masses with the experimental ones. 

The main result of this paper is a new measurement of the quark masses which 
we believe to be more reliable than previous ones. 

The paper is organized as follows. In section 2 we discuss the theoretical 
definitions of lattice quark masses and the renormalization procedures. In 
section 3 we relate lattice quantities to the continuum ones and we review 
the NP renormalization of bilinear operators. In section 4 we give the details 
of the lattices used to extract masses and matrix elements. In section 5 we 
report our results on quark masses in the MS-scheme and finally we give our 
conclusions. 



2 Quark Masses 

The usual on-shell mass definition cannot be used for quarks since they do 
not appear as physical states. Thus the values of the quark masses depend on 
the definition adopted. In the following we will give results for quark masses 
in the MS scheme. 

The quark mass can be defined by the perturbative expansion of the quark 
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propagator renormalized at a scale fi. This is equivalent to the definition from 
the renormalized Vector and Axial Ward Identities, which can be used to give 
a fully non-perturbative determination of the quark mass. 



2. 1 Quark masses from the Vector Ward Identity 



Quark masses can be defined from the Vector Ward Identity (VWI). We first 
consider the Wilson formulation and then describe the Clover Improved case. 
The Vector Ward Identity between on-shell hadronic states a and f3, neglecting 
terms of 0(a), can be written as [15] 



(a\d»V°\f3) = <o# 



A^ 

y 



M 



where ip is a flavour triplet of bare quark fields M is the quark bare mass 
matrix, A a /2 are the generators of the SU(3) flavour group normalized by 
Tr[X a \ b ] = 2 5 ab and the derivative <9 M is the usual asymmetric lattice deriva- 
tive. The conserved vector current V,? is: 



v,?(x) = - 

j 2 



A a 



A a 



+ ^(x + /i) (7^ + 1)^(^) ^(3;) 



(2) 



and its renormalization constant satisfies Zy — 1. are the link variables on 
the lattice. By taking the linear combination = — iV^ we obtain: 



<^W> = i(l 



>a\S\P) 



(3) 



where fcj is the hopping parameter of i-ih quark, S = q 2 qi is the bare scalar 
density. 

On the lattice the ratio of the vector and scalar matrix elements could be used 
to determine quark masses. Unfortunately the scalar matrix element turns out 
to be extremely noisy, preventing any reliable analysis. We then use eq.(3) only 
to fix the relation between the lattice bare quark mass in lattice units and the 
hopping parameter, i.e. 



where k c is the critical value of the hopping parameter. Eq. (3) shows that 
the quark mass, defined in (4), is renormalized with Zg 1 = Z m , where Zg is 
the renormalization constant of the scalar density. For the pseudoscalar and 
vector meson masses as a function of quark masses, ChPT gives 

M 2 PS = CVi + m i) , (5) 
M v = A + B{m l +m 2 ) . (6) 



From equations (5) and (6) one can extract the bare lattice quark mass of a 
given flavour. This is the so called "spectroscopy method" . 



In the Improved case eq. (3) will be more complicate as other terms will appear 
due to our use of "-^-rotated operators" [14] : 



jl - i zlj) - (1 - z)mi J qi 



(7) 



where Ij) is the lattice covariant derivative, i denotes the flavour of the quark 
field and z is an arbitrary real number. For z = 0, each field qi is multiplied 
by a numerical factor so that the tree-level improved bilinear operators are 
simply [16] 

Ol = (l + m t )O r , (8) 



obtaining, as for eq. (3), 



H^|/?> = (9) 
^(mi - m[) (a\{2S\x) + S\x + 6) + S\x - 0)|/3) , 

with the symmetric derivative d^f(x) = (f(x + fi) — f(x — fi))/2a [17], where 
a is the lattice spacing. 

However, one should note that the combination of scalar densities in the right 
hand side of eq. (9) can be replaced by S r (x) up to higher orders in the lattice 
spacing a. From eqs. (8)- (9) and the fact that the Clover Improved action has 
only 0(ag 2 ) errors, we obtain that the mass mf, improved up to 0(ag 2 ), is 

r 1/1 lu 1/1 lxn 

= 2(5 - it;)! 1 -ik -*;)]■ (10) 

Eq. (10) is, of course, the tree level approximation of the relation obtained in 
[18]. This procedure can be straightforwardly applied to the perturbative case 
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using RC's that do not contain the field "rotation" contribution, but not in 
the non perturbative one as, for historical reasons, we have computed RC's 
using "rotated" fields with z — 1. 



2.2 Quark masses from the Axial Ward Identity 



In this section we will describe the definition of the mass from the Axial Ward 
Identity (AWI). As for the VWI case we first consider the Wilson formulation 
and then describe the Clover Improved case. Close to the chiral limit and 
neglecting terms of 0(a), the Axial Ward Identity is [15]: 

Z A (a\d»A«\P) = (o#{y , M - M}l,m , (11) 



where is the RC of the axial current and M, defined in [15], in the chiral 
limit is the critical value of the quark mass. Following the same steps as for 
the VWI, and for the same flavour choice as in eq. (3), we obtain that one can 
extract quark masses from the ratio of matrix elements 

2 {pl+p2) ~ (0\P\PS(p = 0)) ' (12) 



where A^ = o 2 7|u75<Zi, P = Q2I5Q1 and PS is a pseudoscalar meson. This ratio 
is renormalized with Za/Zp, where Zp is the RC of the pseudoscalar density. 

For the Clover improved case, the only difference is that we need to consider 
the ratio of improved operators and to take the symmetric derivative, as in 
eq. (9). 



3 Renormalization of quark masses 



In the continuum, "physical" quantities like for example quark masses and de- 
cay constants are defined specifying the renormalization prescription used to 
eliminate the divergences in the matrix elements of the operators determining 
them. The standard practice consists in choosing a Dimensional Regulariza- 
tion scheme, for instance, NDR — MS. On the lattice, however, the matrix 
elements of the corresponding operators are regularized by introducing a hard 
cut-off, a -1 , in the momentum integrals. Since the continuum and lattice 
renormalization schemes, in general, differ, so do the corresponding ampli- 
tudes. Therefore, in order to convert the bare lattice values of amplitudes to 
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their continuum counterparts, we have to know the matching coefficient, the 
so-called lattice renormalization constant. 



One method, which is described in section 3.1, consists in computing two 
amplitudes to a given order in perturbation theory in the continuum and on 
the lattice. Imposing that they must be equal after renormalization at a given 
scale, one finds the lattice RC's to some order in perturbation theory. 

However, to avoid the use of perturbation theory on the lattice, a general 
method for the non-perturbative calculation of lattice RC's in the RI (regu- 
larization independent) scheme has been devised. We shall briefly discuss this 
in section 3.2. Once we have obtained the lattice RC's in the RI scheme non- 
perturbatively, the conversion to other continuum schemes like NDR — MS 
and the running of the continuum quantity to other scales can be calculated 
in perturbation theory at a given order. 



3. 1 The perturbative approach 



In the perturbative approach one uses lattice and continuum perturbation 
theory at the next-to- leading order (NLO). 

For the VWI definition, we have 



m MS (ri = U^^K/a)Jpma- 



a 8 (n/a) 
1 + 47r Kvwi 



mar v (13) 



where m is the bare quark mass of a generic flavour, 

Kvwi = ~C F (A S + A s ) - 7 (°> logfr) , (14) 



Cf = {N 2 — 1)/2N and N being the number of colours. As and As have been 
computed in ref. [16,19] and following them can be extracted from table 1 for 
both actions. 



/ a 8 (fj) \7 (u) /2ft, 



1 + 



aM-a a (Tr/a) n W 7 (0) A 



47T 



2/3b 2/3 2 



(15) 



where the coupling constant is defined in the MS scheme, 
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UN-2n f 
A)- 3 

n 34 ,r2 10 Ar iV 2 — 1 

A^-^n,-— n, 

7 (0) =6^ (16) 

m 97iViV 2 -l /iV 2 -K2 10n f iV 2 -l 

7 W = + 3( J ^ 

[ 3 2JV MJV ^ 3 2iV 

(17) 

and n/ is the flavour number. In the improved case, as discussed before, As in 
eq. (14) should not include the contribution due to the rotation of the fields. 

Analogously, for the AWI method, we have 



r?MS 

m MS (») = U™ s Ma)%pa-i 



J p 

tMS 



U^Qm,*/ a) 



a s {ix/a) 
J- H : J^awi 

4"7T 



pa' 1 , (18) 



where | is the ratio defined in (12) for a generic flavour and 

K AWI = C F (A A - A P ) - 7 (°) log(Tr) . (19) 

Also A a and A P can be computed from table 1 following [16,19]. 
3.2 The non-perturbative approach 



At scales a -1 ~ 2—4 GeV, where a is the lattice spacing, of our simulations, we 
expect small NP effects on the renormalization constants of bilinear operators. 
However "tadpole" diagrams [20], which are present in lattice perturbation 
theory, can give rise to large corrections and then to large uncertainties in the 
matching procedure at values of (3 — 6/gl — 6.0 — 6.4. These problems are 
avoided using NP renormalization techniques [21,22]. Although the nature of 
the systematics changes in adopting a NP approach, mainly because of 0(a) 
effects as we shall discuss in [23], we believe that the NP determination of 
RC's is more reliable. 



We shortly review the NP method for the RC's which enter the determination 
of the quark masses. For the full discussion of the method, the results obtained 
and their systematics we refer to a forthcoming paper [23]. 
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Let us consider a quark bilinear Or = qFq, where T is a Dirac matrix and q 
now represents the generic lattice flavour-degenerate quark field. In this work 
we will consider the scalar and pseudoscalar densities and the Axial Vector 
current. 

The renormalization conditions are imposed on the amputated Green functions 
computed between off-shell quark states of momentum p in the Landau gauge 

Koipa) = S q {pa)~ l Go{pa)S q {pa)- 1 (20) 



where Goipo) and S q (pa) are the non- amputated Green functions and the 
quark propagator, calculated non-perturbatively via Monte Carlo simulations 
[21]. Possible effects from Gribov copies or spurious solutions [24] have not 
been considered. The RC ZQ l (pa,g ) of Or is determined by the condition 

Z^HVW^ P Ao(pa)\ P *=^ = 1, (21) 



where Po is a suitable projector on the tree- level amputated Green function 
(normalized to 1) [21] and Z q is the wave function RC which can be defined 
in different ways [21]. From the Ward Identities 

(22) 

p 2 =fi 2 



Z q (pa) = -i-Tr( ) 



To avoid derivatives with respect to a discrete variable, we have used 



Z'Ua) 



12 4£^= M sin 2 (^a) 



(23) 



which, in the Landau Gauge, differs from Z q by a finite term of order a 2 s which 
is irrelevant at the NLO. 

For this procedure to be reliable \x must satisfy the condition p, <^ 1/a to avoid 
discretization errors but also p ^> Aq C d to avoid non-perturbative effects or 
higher order corrections in the continuum perturbative expansion. 

The quark mass in the MS scheme is then defined as 

m R V) , (24) 



™ ms M = C(m') 



i + 



a s(n') nLA N 
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Regularization 


A s 


A Loc 


A® 


A Loc 


A® 

L±p 


A Loc 


A® 


Wilson 


-12.85 


-0.10 




-9.78 




-3.0 




Clover 


-9.2 


-10.07 


12.0 


-13.1 


-4.3 


-4.6 


11.7 



Table 1 

Coefficients entering the perturbative matching of the quark masses, 
where [21] 

^lan _ a N 2 — 1 
2A/ 



(25) 



and for m RI (/j,') we get 



m m (ji') 
m RI (fx') 



-ma 



7 RI n 

Zf{ l i>a)T 



(26) 



for the Vector and Axial Vector WI respectively. As both the RI and the MS 
respect chirality we have Z RI = Z/ 



'A- 



This discussion applies immediately to the unquenched case [25]. 



4 Details of the analysis 



In this section we describe the extraction of meson masses, matrix elements 
and lattice quark masses from the two point correlation functions and the 
interpolation/extrapolation of the results in the heavy and light quark masses 
to the physical points. We also try to estimate the overall systematic errors 
on the lattice quark masses. 



4-1 Lattice details 



In this work we have used various lattices that have been generated by the 
APE group in the last years. Tables 2 and 3 show the parameters of the lattices 
that we have analyzed. A more detailed discussion of lattice calibration and 
spectroscopy can be found in [7] for the light quark systems and in [26] for 
the heavy-light ones. Here we will just summarize the most important points. 
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Meson masses and axial-pseudoscalar matrix elements have been extracted 
from two-point correlation functions in the standard way from the following 
propagators 

X 

G 05 (t)=J2(A (x,t)P\0,0)} (27) 

X 

and 

Gu(t)= EEWM)V?(0,0)>, (28) 

i=l,3 x 

where 

P(x,t) =q(x,t)-f 5 q(x,t) , 
A^(x,t)=q(x,t)-f^ 5 q(x,t) , 
Vi(x,t) =q(x,t)7iq(x,t) . 

where q represents the generic lattice flavour-degenerate quark field and % is 
a spatial index. We fit the zero-momentum correlation functions in eqs. (27) 
and (28) to a single particle propagator 



^55 -y 

G 55 (t) = —— exp(--M PS T) cosh[M P5 (- - t) 



G u (t) 



M PS 
z u 



Mi 



v 







) cosh 


M v (^-t) 



(29) 



in the time intervals reported in tables 2 and 3. In (29), T represents the 
lattice time extension, the subscripts PS and V stand for pseudoscalar and 
vector mesons. To improve stability, the meson (axial-pseudoscalar) correla- 
tion functions have been symmetrized (anti-symmetrized) around t = T/2. 
The time intervals for the fits are chosen with the following criteria: we fix 
the lower limit of the interval as the one at which there is a stabilization of 
the effective mass, and, as the upper limit, the furthest possible point before 
the error overwhelms the signal. The errors have been estimated by a jack- 
nife procedure, blocking the data in groups of 10 configurations and we have 
checked that there are no relevant changes in the error estimate by blocking 
groups of configurations of different size. 



10 





C60a 


C60b C60c C60d 


C62a 


C64 


a 
P 


6.0 


6.0 6.0 6.0 


6.2 


6.4 


jp L-onts 


490 


600 200 200 


250 


400 


Volume 


18 3 x 64 24 3 x 40 18 3 x 32 16 3 x 32 24 3 x 64 24 3 x 64 


Kl 


0.1425 


0.1425 0.1425 0.1425 


0.14144 


0.1400 




0.1432 


0.1432 0.1432 0.1432 


0.14184 


0.1403 




0.1440 


0.1440 0.1440 0.1440 


0.14224 


0.1406 








0.14264 


0.1409 


u 






0.1210 










0.1250 










0.1290 










0.1330 








Light-light mesons 






h - t 2 


15-28 


15-19 11-15 11-15 


18-28 


24-30 






Heavy-light mesons 






h-h 






20-28 




a-\K*) 


2.12(6) 


2.16(4) 2.07(6) 2.23(9) 


2.7(1) 


4.0(2) 



Table 2 

Summary of the parameters of the runs with the SW-Clover fermion action analyzed 
in this work. 



We extract p from the ratio 



P = 



( d A (t)P(0)) 
(P(t)P(O)) 



(30) 



at zero spatial momentum. At large time separations, we fit the ratio G^{t) /G^{t) 
and we get 



p = sinh(Mp S ') 



(Mt)P(Q)) 

(P(t)P(O)) 



(31) 



in a t region where the signal stabilizes. The hyperbolic sine comes from the 
discrete symmetric derivative which is necessary to consistently take into ac- 
count O(a) terms. 



Tables 2 and 3 also contain the values of the lattice spacing a extracted from 
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WbU Wbza Wdzd Wo4 




# Confs 
Volume 


6.0 6.2 6.2 6.4 
320 250 110 400 

18 3 x 64 24 3 x 64 24 3 x 64 24 3 x 64 


h 


0.1530 0.1510 0.1510 0.1488 
0.1540 0.1515 - 0.1492 
0.1550 0.1520 0.1520 0.1496 
0.1526 0.1526 0.1500 


hi 


0.1255 0.1300 0.1300 
0.1320 0.1350 0.1350 

0.1385 0.1400 0.1400 

n i /ion n 1 /i n 1 a 
U.14ZU U.140U U.140U 

0.1455 0.1500 


u - u 


Light-light mesons 
15-28 18-28 18-28 24-30 


h-t 2 


Heavy-light mesons 
15-28 20-28 20-28 


a- l {K*) 


2.26(5) 3.00(9) 3.0(1) 4.1(2) 



Table 3 

Summary of the parameters of the runs with the Wilson action analyzed in this 
work. 

the K* mass with the "lattice physical plane" (lp-plane) method [7]. 

We have also performed dedicated runs at f3 — 6.0, 6.2 and 6.4 for both the 
Wilson and the SW-Clover actions to measure the RC's according to the 
discussion of section 3.2. The parameters of these runs and the results for 
Zg 1 , Zp 1 and Z^ 1 are reported in table 4. 

4-2 Extraction of raw results for light and strange quarks 

Once the hadronic correlation functions have been fitted, and the lattice 
masses and matrix elements extracted, we have to perform a number of inter- 
polations/extrapolations to extract physical quantities. 

We extract light and strange quark masses from the meson spectroscopy and 
from the Axial Ward Identity with three different methods: 
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n 
O.U 


r n 
O.U 


D.Z 


D.Z 


0.4 


ft A 


Action 




o W 


Wilson 


O W 


Wilson 


Q.XKT 
O W 


Wilson 


# Confs 




1 nn 
1UU 


1 nn 
1UU 


loU 


1 nn 
1UU 


DU 


ftn 
DU 


Volume 




16 3 x 32 


16 3 x 32 


16 3 x 32 


16 3 x 32 


24 3 x 32 


24 3 x 32 


k 




(J. 1425 


n -i r on 


n A A A A A 

U. 14144 


(J.lolU 


n i a nn 

(J.14UU 


n 1 /ioo 








n 1 ^/in 

U. lO^U 




U. 1010 


U. 14iUo 








0.1440 


0.1550 


0.14224 
0.14264 


0.1520 
0.1526 


0.1406 
0.1409 


0.1496 
0.1500 


k c 


0.14551 


0.15683 


0.14319 


0.15337 


0.14143 


0.15058 


Zg 1 \m q a 


= 0) 


0.834(18) 


0.682(9) 


0.851(11) 


0.722(8) 


0.852(13) 


0.742(8) 


ZfP(m q a 


= 0) 


0.409(8) 


0.447(5) 


0.466(4) 


0.499(5) 


0.555(6) 


0.572(4) 


Zf[m q a 


= 0) 


1.047(18) 


0.808(7) 


1.023(4) 


0.812(6) 


1.012(9) 


0.825(6) 



Table 4 

Parameters of the runs used in the non perturbative calculation of the renormal- 
ization constants and values of Zg 1 , Zp 1 and Z^ 1 at a scale \ia ~ 1. The errors 
reported are statistical only. 

(1) We define the physical plane [My, M| s ] (lp-plane [7]) and, assuming 
that only linear terms are important, we determine the lattice values 
and M K * imposing that M V /M PS coincides with the experimental value 
C s i = Mk?/Mk P [27]. Comparing M K *a~ x with its experimental value 
we extract the lattice spacing reported in tables 2 and 3. Using relations 
(5) and (4) and imposing the prediction of the CPTh [1] 



rji 

R= = = 24.4 ±1.5 , 

m 



(32) 



where m = (m u + m,d)/2, we obtain 



m. 



m - 



C(l + l/R) 
R ' 



(33) 
(34) 



The lattice strange quark mass obtained with (33) for all lattices are 
reported in table 6. 

In extracting the values of a -1 , M K and C from our data, we assume 
that the relation between quark masses and squared pseudoscalar masses 
is linear up to the strange quark region (see eq.(2)), so we can compute 
these quantities using our data for quark masses in the "strange mass 
region" by interpolating only. In other words, no extrapolation to the light 
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quark region is performed. Possible deviations from linearity induced by 
terms of the form D (m s + m) 2 + D' (m s — m) 2 are known to be small 
and will give rise to uncertainties which can be neglected. Moreover, the 
value of k c is not required to obtain the quark masses in the Wilson case. 
For the improved action, however, a small dependence of C on k c exists 
because of the quadratic term in eq.(7). This gives rise anyway to quite 
small fluctuations. 

For the Axial Ward Identity, in the same fashion, we define the plane 
[p, M PS ], where p is defined in the equation (31), and we obtain 

P\m ps =M k /,,\ 

^ = (TTW (35) 

P=f- (36) 

In this case, for both actions, there is no chiral extrapolation needed. The 
values of p s for all lattices are reported in table 7. 
(2) On the physical plane [My Mp s ], we determine M K and M K * as in (1) 
and M n and M p imposing that My/Mps coincides with the experimental 
value Cu = M e p xp /M* xp . Using (5) and (4) we obtain 

M 2 

2 < 37 > 

m. = i(<-^). (38) 

The lattice light quark masses obtained with (37) for all lattices are 
reported in table 6. For the Axial Ward Identity we define the plane 
[p, Mp S ] and we obtain 



1 

P = ~P 

^ M PS =M lr 

1 

Ps — P\Mps=M k ~ ^ 



(39) 

M PS =M 1T 



The values of p for all lattices are reported in table 7. 
(3) Using eqs. (4), (6) and the lattice spacing determined as in (1) we compare 
(M^a^ 1 ) to its experimental value and we obtain 

Mi xp a-A 

m s = . (40) 

From (4) and (40) we determine k s and then we use 

9 = G(i - y) (41) 



to determine p s 
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The three methods (1), (2) and (3) described above should give consistent 
values for m, m s , p and p s , apart from discretization errors and quenching 
effects. It is well known [7,10] that the strange quark mass obtained from 
(3) is systematically higher than the value obtained from M K (1). 

The ratio of the strange quark masses obtained with different methods is 
connected to experimental quantities. From (5) and (6) we get 

M v = A + DM 2 PS (42) 



where D ■ C = B. A and D can be measured fitting linearly the experimental 
data on the meson masses [27]. To compare lattice data with experimental 
measurements we prefer to use two equivalent dimensionless quantities 



< 43 > 

In table 5 we report the experimental values and the lattice predictions we 
have obtained for J and L. One can see that J has quite large fluctuations 
and is also quite far from the experimental value. It has been argued [28] that 
this may be mainly due to the quenching approximation. 



Run 


J 


L 


Rm s 


Exp. 


0.499 


0.847 


1 


C60a 


0.367(13) 


0.887(4) 


0.852(17) 


C60b 


0.377(8) 


0.884(2) 


0.866(10) 


C60c 


0.359(16) 


0.890(5) 


0.842(22) 


C60d 


0.393(11) 


0.879(3) 


0.884(14) 


W60 


0.339(10) 


0.896(3) 


0.815(14) 


C62a 


0.343(50) 


0.895(15) 


0.821(68) 


W62a 


0.360(21) 


0.889(7) 


0.843(28) 


W62b 


0.380(28) 


0.883(9) 


0.869(35) 


W64 


0.386(16) 


0.881(5) 


0.876(19) 


C64 


0.401(17) 


0.877(5) 


0.894(20) 



Table 5 

Experimental and lattice values for J, L and R ms . 
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The ratio of the quark masses obtained with different methods is only function 
of J and L and experimental ratios. In particular the ratio of the strange quark 
mass obtained with (1) and (3) is 

R = (!^i£ = I 2 J__ m 



where R is defined in eq. (32). In Table 5 we report the experimental value 
for R^ p and the values we have obtained for our lattices. R ms is a ratio of 
different definitions of lattice quark masses that does not suffer from 0(a) 
effects induced by renormalization constants and is related to experimental 
quantities. We believe that the difference R^ xp — R^ * s a good estimate of 
the overall systematic error on the strange quark mass. This error should 
essentially take into account both the quenching approximation and the 0(a) 
effects that we cannot a priori estimate. 



4-3 Extraction of raw results for the charm quark mass 



In order to obtain the charm quark mass we have to extrapolate the meson 
masses and the matrix elements both in the heavy and light quark masses. It 
is clear that in this case 0(ma) effects will be much more relevant than for 
the strange and light quarks. We use two different methods: 

(1) We first interpolate/extrapolate linearly the heavy-light pseudoscalar me- 
son masses in the light quark mass to m s extracted from (33). We then 
determine the charm quark mass by fixing the D s -meson mass to its phys- 
ical value and using the equation (4). The lattice charm quark masses 
obtained in this way for all lattices are reported in table 6. For the Axial 
Ward Identity p c is given by 

Pc = p\m ps =m Ds - Ps (45) 

with p s defined in eq. (35). The values obtained for all lattices are reported 
in table 7. 

(2) As in (1) but we determine the value of the charm quark mass by fixing 
the D*-meson mass to its physical value. 

As for the light quark masses, the two methods described above should give 
consistent results for m c and p c , apart from discretization errors and quenching 
effects. We shall use the spread between the two determinations as an estimate 
of the overall systematic error. 
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5 Physical Results for quark masses 

On the basis of the discussion of previous sections, we now present the final 
results. From the lattice data of the left columns of tables 6 and 7, we obtain 
the PT and NP MS results reported in the same tables. 

In the PT case there is a large ambiguity in the choice among the lattice strong 
coupling constants which have been proposed to have a better convergence of 
the series. In the appendix we give in detail which definitions we have used. The 
difference among them is of course 0{a 2 s ) but on the lattice the behaviour of 
the perturbative series can be very bad [20]. The different values are averaged 
and the spread taken as an estimate of the error. 

In the continuum we have used both in the PT and NP case eq. (47) with 
iif = and Aqcd — 251 ±21 MeV [29] as discussed in the appendix. We have 
checked that the final results change by less then 3% if we had used rif — 4 
and Aqcd = 340 ± 120 MeV, and this is only part of the error due to the 
quenching approximation which we cannot estimate. 

The NP results contain only the statistical error on the matrix elements and 
renormalization constants. The AWI results are shown in figs. 1-4. 

In the spectroscopy case we find a reasonable agreement between the NP 
and PT results that turn out to be compatible with previous determinations 
[6,7,10]. This is not the case for the AWI where the PT results are lower 
than the NP ones by more than two standard deviations. In our opinion this 
confirms that the perturbation theory fails in the determination of the pseu- 
doscalar RC. On the other hand the NP results for different actions and meth- 
ods are in very good agreement among them. 

As far as the a dependence is concerned, for (3 = 6.0 and 6.2 the data are 
definitely stable, within the errors, for all methods. At (3 — 6.4 the data 
apparently show a slight decrease. Due to the small physical volume at this (3 
value and the fact that this effect is present both in the Wilson and the Clover 
case we believe that one cannot disentangle volume and 0(a) effects in this 
lattice [7] and we will not take into account in any further result. With our 
present data an extrapolation in a is then out of reach. 

The charm quark results appear to be quite more noisy than the light and 
strange ones. This supports the fact that at large masses order am contami- 
nations are important. 

To obtain final results, we average the non perturbative AWI results at /3 — 6.0 
and 6.2 as independent ones. For the strange and charm quark, we also take 
into account the overall systematic error which we evaluate from the spread in 
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the quark masses extracted from different mesons as described in sec. 4. This 
error is also propagated to the light quark using eq. (34) and is the second 
one in eq. (46). 

From NP renormalized AWI we obtain our final results 
m^(2 GeV) = 5.7 ± 0.1 ± 0.8 Mev 

mf^(2 GeV) = 130 ± 2 ± 18 Mev (46) 
mf*(2 GeV) = 1662 ± 30 ± 230 Mev . 



which are reported in the abstract. 



Run 


ma 1 


m' 


MS 


m s a 1 


m 


MS 
s 


m c a 1 


mf s 






NP 


Pert. 




NP 


Pert. 




NP Pert. 


C60a 


3.45(12) 




4.9(4) 


78.9(22) 




113(8) 






C60b 


3.36(9) 




4.8(3) 


77.3(17) 




110(7) 






C60c 


3.49(14) 




5.0(4) 


79.4(25) 




113(8) 






C60d 


3.21(13) 




4.6(3) 


74.7(26) 




107(7) 






W60 


4.36(12) 


5.8(2) 


6.1(8) 


97.9(20) 


129(4) 


137(17) 


1335(31) 1764(62) 1865(247) 


C62a 


3.60(29) 




5.2(6) 


81.1(41) 




117(9) 


881(18) 


1275(90) 


W62a 


4.10(18) 


5.6(3) 


5.9(8) 


93.3(29) 


127(5) 


134(17) 


1205(24) 1646(44) 1730(212) 


W62b 


3.98(23) 


5.5(4) 


5.7(8) 


91.8(39) 


126(6) 


132(17) 


1206(32) 1651(64) 1738(222) 


W64 


3.55(18) 


5.1(3) 


5.3(7) 


82.0(35) 


119(6) 


123(17) 






C64 


2.86(17) 




4.3(4) 


66.9(34) 




100(9) 







Table 6 

Lattice quark masses and the corresponding MS values in MeV at NLO from the 
pseudoscalar meson spectroscopy and a -1 from Mk*. MS masses are at a scale 
H = 2 GeV. 



6 Conclusions 

We have discussed the quark mass renormalization. We have calculated the 
quark masses from the meson spectroscopy and from the Axial Ward Identity 
using different sets of quenched data with /3 = 6.0, 6.2 and 6.4 and using 
the Wilson and the "improved" SW-Clover action. The data at (3 = 6.4 have 
been used only for an exploratory study as the physical volume and the time 
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Run 


pa 1 


m 


MS 


p s a 1 


m 


MS 








fN.r 


.ren. 






r^ert. 


IN r r^ert . 


LbOa 


2.b4(ll) 


b.0(3) 


/i n / o \ 

4.0(3) 


b0.4(20) 


13b(b) 


93(8) 




CbOb 


2.52(7) 


5.7(2) 


o r\ / o \ 

3.9(3j 


58.0(14) 


131(5) 


on / \ 

89(b) 




CbUc 


2.b8(ll) 


b.0(3) 


/i 1 /o\ 

4-1(3) 


b0.9(20) 


136(6) 


93(7) 




C60d 


2.41(12) 


5.5(o) 


3.7(3j 


55.1(24) 


128(5) 


86(7) 




WbO 


Q /I O /l 1 \ 

o.4o(ll) 


0. 7 (2) 


4.3(5 J 


f 5.2(21) ) 


12 / (5) 


no / 1 1 \ 

yo(ll) 


1091(30) 1777(07) 13b2(lb0) 


Lb2a 




b.0(5) 


4.4(5j 


b3.b(33) 


10/1/ o\ 

134(8) 


99(9) 


820(3b) 1725(93) 12b9(113) 


VV OZa 








78 Qf9"^ 


1 9fif ^ 

IZD^O J 






W62b 


3.32(19) 


5.3(3) 


4.3(5) 


76.5(32) 


124(6) 


99(12) 


1001(29) Ibl0(b6) 1301(149) 


W64 


3.19(17) 


4.9(3) 


4.3(5) 


73.7(33) 


114(6) 


100(13) 




C64 


2.41(14) 


4.7(3) 


3.8(4) 


56.3(27) 


110(6) 


90(8) 





Table 7 

Lattice quark masses and the corresponding MS values at NLO, in MeV from the 
Axial Ward Identity and a -1 from Mk*- MS masses are at a scale p = 2 GeV. 



extension of the lattice may be too small to be reliable. 

Perturbation theory appears to fail giving inconsistent results for the two 
methods. The results for independent NP methods are well consistent and 
stable strongly supporting the reliability of our final results reported in the 
abstract. In the f3 range that we have considered and with our statistics, we 
do not believe that one can safely extrapolate to the continuum limit both for 
the Wilson and the SW-Clover action. 



7 Appendix 



In this appendix we shall describe the different definition of the lattice strong 
coupling constant that we have used for the PT results. 



The running coupling constant in MS scheme is: 



«4 2 ) = 1 
4tt PM» 2 /^qcd) 



+ 



(47) 
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this is related to the bare lattice coupling (a) = g 2 L (a)/An by the relation: 



o£(a) 



a s (q) 



1 + 



An 



7T 



A,ln(— ) 2 + 48.76 
ga 



(48) 



In ref. [20] a new expansion parameter, a^ 1 has been introduced based on the 
heavy quark potential: 



o£(a) 



ct 



Vi 



' 1+ «- Vl (9) 



47T 



7T 



/3 ln(—) 2 + 59.09 
ga 



(49) 



using the perturbative expansion for the trace of the plaquette one obtains, 
[20]: 



(a) 



a 



a s 

1 + — 



ir /a) 



An 



6.45 



(50) 



where Up is the plaquette. The third definition that we have used is: 

i my 



a 



ia) 



a 



v 2 



(51) 



where k c is the critical value of the hopping parameter. 



In our analysis we have used a s , a^ 1 , a 2 as they have been used in the 
literature in the past. One should note that each definition has different 0(a) 
effects and eq. (51) is only valid at the LO while eqs. (47-50) are valid up to 
NLO. We have made this choice as it is the same done in ref. [6] and so the 
results can immediately be compared. 

We have evaluated all cases both in the quenched and unquenched case. In 
the latter case the value of A.Q f CD at the various scales is obtained from the 
SLD-LEP result A b QCD = 240 ± 90 MeV [27] by matching the strong coupling 
constant at n/ and rif + 1 at the threshold of the new flavour [30]. In the 
quenched case we have used Aqcd = 251 ± 21 MeV [29]. Furthermore we 
have also chosen two different scales used to evaluate the strong coupling 
constant, i.e. a s {it/a) and a s (l/a). All these results have been averaged and 
the spread taken as the error. We have also done this analysis using Aq C d = 
340 ± 120 MeV and we find that the final average changes by less then 5% 
which is well within the error. As an example example af s changes from 0.20 
to 0.31 at a scale fj, — l/a — 2 GeV which leads to a change in the correction 
of that order. 



20 



Acknowledgements 

We wish to thank the APE collaboration for allowing us to use the lattice cor- 
relation functions presented here. We warmly thank C. R. Allton, V. Lubicz, 
E. Franco, G. Martinelli, M. Testa and A. Vladikas for enlightening discus- 
sions. We also acknowledge the use of C. R. Allton's analysis program. M.T. 
acknowledges the support of PPARC through grant GR/L22744.C 

References 

[I] H. Leutwyler, Nucl. Phys. B 337, (1990) 108 and references therein. 
H. Leutwyler, Phys. Lett. B 378, (1996) 313. 

[2] J. Bijnens, J. Prades and E. de Rafael, Phys. Lett. B 348 (1995) 226. 

[3] M. Jamin and M. Miinz, Z. Phys. C 66 (1995) 633. 

[4] P. Colangelo, F. De Fazio, G. Nardulli, N. Paver, Phys. Lett. B 408 (1997) 340. 

[5] C.A. Dominguez, Proceedings of the 16th International Workshop on Weak 
Interactions and Neutrinos (WIN 97), Capri, Italy, 1997. 

[6] C. R. Allton et al., Nucl. Phys. B 431 (1994) 667. 

[7] C.R. Allton, V. Gimenez, L. Giusti and F. Rapuano, 
Nucl. Phys. B 489, (1997) 427. 

[8] JLQCD Collaboration (S. Aoki et al.), Nucl. Phys. B (Proc. Suppl.) 63 (1998) 
275. 

[9] CP-PACS Collaboration (S. Aoki et al.), 
Nucl. Phys. B (Proc. Suppl.) 63 (1998) 161, 
Nucl. Phys. A (Proc. Suppl.) 60 (1998) 14. 

[10] T. Bhattacharya and R. Gupta, Phys. Rev D 55 (1997) 7203, 
Nucl. Phys. B (Proc. Suppl.) 63 (1998) 95. 

[II] B. J. Gough et al., Phys. Rev. Lett. 79 (1997) 1622. 

[12] SESAM-Collaboration (N. Eicker at al.), Phys. Lett. B 407 (1997) 290. 

[13] B. Shcikfiolcslami, R. Wohlert, Nucl. Phys. B259 (1985) 572. 

[14] G. Heatlie et al., Nucl. Phys. B352 (1991) 266, Nucl. Phys. B (Proc. Suppl.) 17 
(1990) 607. 

[15] L.H. Karsten and J. Smit, Nucl. Phys. B 183 (1981) 103. 
M. Bochicchio et al., Nucl. Phys. B 262 (1985) 331. 



21 



[16] E. Gabrielli et al, Nucl. Phys. B 362 (1991) 475. 

A. Borrelli, C. Pittori, R. Frezzotti and E. Gabrielli, Nucl. Phys. B 409 
(1993) 382. 



[17! 

[is; 

[19 

[20! 
[21 
[22 

[23! 

[24 

[25 

[26! 
[27! 
[28! 
[29! 
[30 



M. Crisafulli, V. Lubicz and A. Vladikas, Eur. Phys. J. C 4 (1998) 145. 

M. Liischer, S. Sint, R. Sommer and P. Weisz, Nucl. Phys. B 478 (1996) 365. 

G. Martinelli and Y.C. Zhang, Phys. Lett. B 123 (1983) 433. 

H. W. Hamber and C. M. Wu, Phys. Lett. B 133 (1983) 351. 

G.P. Lepage and PB. Makenzie, Phys. Rev D 48 (1993) 2250. 

G. Martinelli et al, Nucl. Phys. B 445 (1995) 81. 

K. Jansen et al., Phys. Lett. B372 (1996) 275. 
M. Liischer et al., |hep-ph/9711205| . 

V. Gimenez, L. Giusti, F. Rapuano, M. Talevi, hep-lat /9806006 , to be published 
in Nucl. Phys. B. 

C. Parrinello, S. Petrarca and A. Vladikas, Phys. Lett. B 268 (1991) 236. 
L. Giusti, Nucl. Phys. B 498 (1997) 331. 

TxL Collaboration (N. Eicker et al.), Nucl. Phys. B (Proc. Suppl.) 60A (1998) 
311 and work in progress. 

C. R. Allton et al., Phys. Lett. B405 (1997) 133. 
Particle Data Group, Phys. Rev D 54 (1996) 1. 

UKQCD Collaboration (P. Lacock et al.), Phys. Rev D 52 (1995) 5213. 
S. Capitani et al., Nucl. Phys. B (Proc. Suppl.) 63 (1998) 153. 
W. Marciano, Phys. Rev D 29 (1984) 580. 



22 



> 



6 — 



5 — 



4 — 



n 1 1 1 1 r 



5] 



<> 



W60 O 
W62a □ 
W62b + 
W64 W 



_l I I I I I I I I I I I I I L_ 



0.2 



0.3 0.4 
a [GeV" 1 ] 



0.5 



0.6 



Fig. 1. Non-perturbatively renormalized quark masses m MS for all Wilson lattices 
from the Axial Ward Identity 
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Fig. 2. Non-perturbatively renormalized quark masses m MS for all SW-Clover lat- 
tices from the Axial Ward Identity 
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Fig. 3. Non-perturbatively renormalized quark masses for all Wilson lattices 

from the Axial Ward Identity 
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Fig. 4. Non-perturbatively renormalized quark masses m^ IS for all SW-Clover lat- 
tices from the Axial Ward Identity 
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